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ABSTRACT 

Smoothed monthly mean coronal mass ejection (CME) parameters (speed, acceleration, cen- 
tral position angle, angular width, mass and kinetic energy) for Cycle 23 are cross-analyzed, 
showing a high correlation between most of them. The CME acceleration (a) is found to be 
highly correlated with the reciprocal of its mass (M), with a correlation coefhcient r = 0.899. 
The force [Ma) to drive a CME is found to be well anti-correlated with the sunspot number (i?z), 
r = —0.750. The relationships between CME parameters and i?z can be well described by an 
integral response model with a decay time scale of about 11 months. The correlation coefficients 
of CME parameters with the reconstructed series based on this model (rfi — 0.886) are higher 
than the linear correlation coefficients of the parameters with i?z (ro = 0.830). If a double decay 
integral response model is used (with two decay time scales of about 6 and 60 months), the 
correlations between CME parameters and i?z improve (rt2 = 0.906). The time delays between 
CME parameters with respect to Rz are also well predicted by this model (19/22 = 86%); the 
average time delays are 19 months for the reconstructed and 22 months for the original time 
series. The model implies that CMEs are related to the accumulation of solar magnetic energy. 
The relationships found can help to understand the mechanisms at work during the solar cycle. 

Subject headings: Coronal Mass Ejections (CMEs); Solar cycle; Sunspots; Time variations, secular and 
long term 



1. Introduction 

A coronal mass ejection (CME) is a striking 
manifestation of solar activity seen in the solar 
corona {e.g., Gosling 1993; Webb & Howard 1994). 
In a typical CME, billions of tons of solar magne- 
tized plasma with energy above 10^^ erg can be 
pushed into the space {e.g., Vourlidas et al. 2000; 
Gopalswamy et al. 2000; Falconer et al. 2002). 
CMEs may cause strong interplanetary distur- 
bances and geomagnetic storms (Gonzalez et al. 
1994; Gopalswamy 2010). Very high-energy (GeV) 
particles generated by CME-driven shocks with 
large Mach numbers are hazardous to our modern 
highly technological equipments {e.g., Roussev et 
al. 2003; Lee 2005). 

CMEs originate from large-scale closed mag- 
netic field structures related to active regions with 
or without filaments or to quiescent filaments 



(Munro et al. 1979; Chen & Shibata 2000; Forbes 
et al. 2006; Gopalswamy 2006). The magnetic field 
in the lower corona is a key element in the genesis 
of CMEs {e.g., Forbes et al. 2006). The fastest 
CMEs originate because of an instability of strong 
magnetic fields (> 100 G), generally in active re- 
gions with sunspots (Falconer et al. 2002; Gopal- 
swamy et al. 2003, 2010) while the high latitude 
CMEs are mainly related to quiescent prominence 
eruptions (Gopalswamy et al. 2003). Sunspots 
represent one of the most obvious manifestations 
of solar magnetic activity (Moradi et al. 2010) 
and its number and characteristics can be con- 
sidered as a measure of the energy supplied to 
the corona (Temmer et al. 2003). Studying the 
relationships between CME parameters and solar 
activity, represented by the international relative 
sunspot number {Rz), can be useful to understand 
the mechanisms at work during the solar activity 
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cycle (Sakurai 1976; Antiochos ct al. 1999; Rust 
2003; Wang et al. 2000; Mittal & Narain 2010; 
Kilcik et al. 2011). 

Since the first CME was discovered in data 
from the Naval Research Laboratory (NRL) coro- 
nagraph mounted on the Orbiting Solar Obser- 
vatory-? (OSO-7) NASA satellite (Koomen et al. 
1975), the properties of CMEs have been care- 
fully examined, particularly after the advent of 
the Large Angle and Spectrometric Coronagraph 
(LASCO, Brueckner et al. 1995; Domingo et al. 
1995) on board the Solar and Heliospheric Obser- 
vatory (SOHO). 

Some CME parameters follow the solar cycle 
rather well while others do not (Gopalswamy et 
al. 2003; Kane 2006; Ivanov et al. 2009; Gopal- 
swamy 2010; Gopalswamy et al. 2010; Gerontidou 
et al. 2010, and references therein). The number 
of CMEs per day shows a strong dependence on 
the solar cycle (Hildner et al. 1976; St. Cyr et al. 
2000; Gopalswamy et al. 2003; Gopalswamy 2006, 
2010; Cremades & St. Cyr 2007; Gerontidou et 
al. 2010). The CME speeds have also been found 
to follow the solar cycle (Gopalswamy et al. 2003; 
Gopalswamy 2010; Cremades & St. Cyr 2007). 
However, the CME mass and angular width do 
not exhibit a significant variation with the solar 
cycle (Cremades & St. Cyr 2007). The number 
of CMEs was shown to follow well the sunspot 
number only during the rising phase of Solar Cy- 
cle 23, while there were large fluctuations in the 
maximum and the declining phase of this solar cy- 
cle (Gerontidou et al. 2010; Gopalswamy 2010). 
The CME sources associated with active regions 
follow the butterfly diagram and appear at lower 
latitudes as the cycle progresses, while those re- 
lated to filaments outside active regions tend to 
migrate towards higher latitudes (Gopalswamy et 
al. 2003; Cremades & St. Cyr 2007; Gopalswamy 
et al. 2010, and references therein). 

It is well known that solar flares (e.g., Aschwan- 
den 1994; Temmer et al. 2003), CMEs (Gonzalez 
& Tsurutani 1987; Tsurutani et al. 2006; Ramesh 
2010; Kilcik et al. 2011), and geomagnetic activi- 
ties (Wilson 1990; Echer et al. 2004; Ramesh 2010; 
Kilcik et al. 2011) often lag behind sunspot ac- 
tivity (i?z) from several months to a few years. 
To have a better understanding of the relation- 
ship between solar (i?z) and geomagnetic activ- 
ity {aa index), Du (2011a) proposed an integral 



response model to describe better this relation- 
ship. It is found that the correlation between 
aa and yielded by this model is much higher 
than that found by a simple point-point corre- 
spondence. Some phenomena can be explained 
by this model, such as the significant increase in 
the aa index over the 20th century (Feynman & 
Crooker 1978; Clilverd et al. 1998; Demetrescu & 
Dobrica 2008; Lukianova et al. 2009), the longer 
lag times at solar (cycle) maximum than at solar 
minimum (Wilson 1990; Wang et al. 2000; Echer 
et al. 2004), and the decreasing trend in the cor- 
relation between aa and Rz (Borello-Filisetti et 
al. 1992; Mussino et al. 1994; Kishcha et al. 1999; 
Echer et al. 2004; Du 2011b). 

This paper investigates in detail the relation- 
ships and time delays between CME parameters 
and sunspot numbers (Rz). The CME parame- 
ters (speed, acceleration, central position angle, 
angular width, mass and kinetic energy) and the 
methods used in this paper are described in Sec- 
tion 2. The linear cross-corrcilations among the 
above parameters and their dependencies on R^ 
are first analyzed in Section 3.1. An integral re- 
sponse model (Du 2011a) is then used to study the 
relationships and time delays between the CME 
parameters and R^ in Section 3.2. In this model, 
the output depends not only on the present input, 
but also on past values. In a second step, this 
model is modified and the relationships between 
the CME parameters and R^ improve (see Section 
3.3). Finally, our conclusions are discussed and 
summarized in Section 4. 

2. Data and Methods 

The data used in this study comprise monthly- 
mean sunspot numbers (Rz) produced by the 
Solar Influences Data Analysis Center (SIDC, 
http:/ /www. sidc.be/sunspot-data/) and CME pa- 
rameters obtained from the SOHO/LASCO CME 
catalog identified since 1996 (http://cdaw.gsfc.- 
nasa.gov/, see e.g., Gopalswamy et al. 2000; 
Yashiro et al. 2004). The following CME parame- 
ters are analyzed: 

1. V\: linear speed obtained by fitting a 
straight line to the height-time measure- 
ments (in units of km s"-'^). 

2. V20- initial speed obtained by fitting a 
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parabola to the height-time measurements. 

3. V2: final speed at the time of final height 
measurements by fitting a parabola. 

4. V20R- speed obtained as above but evaluated 
when the CME is at a height of 2Oi?0 (solar 
radii). 

5. a: acceleration obtained by fitting a parabola 
to the height-time measurements (in km 

6. W: angular width in the plane of the sky (in 
degrees). 

7. Pq: central position angle (CPA, in degrees), 
from solar north (counterclockwise). 

8. Pm- position angle (in degrees) at the fastest 
portion of the leading edge (Yashiro et al. 
2004). 

9. M: mass (in grams). 

10. E: kinetic energy (in ergs), obtained from 
the linear speed (Vi) and the CME mass 
(M). 

The above parameters are first integrated over 
each day and, then, averaged over each month to 
obtain the monthly-means of integrated daily pa- 
rameters. Such a treatment aims to analyze the 
total contribution of CMEs over a month. Ide- 
ally, CME projection effects should be corrected 
in advance to accurately measure the parameters 
(Cremades & St. Cyr 2007; Ivanov et al. 2009; 
Gopalswamy 2010). As doing this for all CMEs 
is a difficult task and there is a small number of 
CMEs with W > 120° (11%, Gopalswamy et al. 
2010), projection effects have not been corrected in 
this study. In addition, we do not consider the N- 
S asymmetry of the position angles (Pq and Pm)', 
that is, we use the values of position angles in 
terms of co-latitude — computed from the solar 
North pole for the northern lieinispliere and from 
the solar South pole for the southern hemisphere. 
The Po and Pm values may be greater than 90° due 
to the integral contributions of many CMEs in a 
day. To filter out high frequency variations in the 
data, the parameters are smoothed with the stan- 
dard 13-month running mean technique (with half 
weights at the two ends). The results are shown in 
Fig. 1 for the period July 1996 to December 2008 
(the time of solar minimum preceding Cycle 24). 



3. Results 

3.1. Correlations between CME Petrame- 

ters 

The relationships among CME parameters and 
other phenomena have been well examined in the 
past. For example, faster CMEs tend to be wider 
{e.g., Gopalswamy 2010; Gopalswamy et al. 2010), 
with a correlation coefficient of r = 0.44 be- 
tween Vi and W for Vi > 900 km s~^ (Yashiro 
et al. 2004), and wider CMEs generally have 
greater mass contents, M oc W^-^ {e.g., Gopal- 
swamy 2010; Aarnio et al. 2011). CMEs with 
the above-average speeds (466 km s~^) deceler- 
ate, while those with speeds well below the aver- 
age accelerate (Gopalswamy et al. 2000; Yashiro 
et al. 2004; Gopalswamy 2006; Gopalswamy et al. 
2010; Aarnio et al. 2011). Strong CMEs are often 
accompanied by solar flares, type II radio bursts, 
energetic particles, and geomagnetic storms {e.g., 
Gonzalez et al. 1994; Roussev et al. 2003; Gopal- 
swamy 2010). 

In this section, we simply analyze the correla- 
tion coefficients between the CME parameters, as 
listed in Table 1. It is seen that these values are 
very high, from 0.77 to 0.99, except for related to 
a, Pq and Pm. The following can be noted: 

1. Vi (solid line), V20 (dotted line), V2 (dashed 
line), and V20R (dash-dotted line) in Fig. la 
are highly cross-correlated with correlation 
coefficients in the range from r = 0.956 to 
0.998. So, if one of these speeds is well cor- 
related with any other parameter {e.g., R^), 
the others will also be. 

2. a (solid line in Fig. lb) is inversely correlated 
with Vx, V20, V2, and V20R, r = —0.436 to 
—0.696, implying that CMEs with faster 
speeds tend to decelerate faster (Gopal- 
swamy et al. 2000; Yashiro et al. 2004; 
Gopalswamy 2006; Gopalswamy et al. 2010; 
Aarnio et al. 2011). 

3. W (solid line in Fig. Ic) is well correlated 
with Vi, V2Q, Vi, and ^zoR, r = 0.820 
to 0.932, suggesting that CMEs with faster 
speeds tend to be wider (Yashiro et al. 2004; 
Gopalswamy 2010; Gopalswamy et al. 2010). 

4. W is well anti-correlated with a, r = —0.805, 
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Fig. 1. — (a) Smoothed monthly moan Vi (sohd hnc), V20 (dotted hne), V2 (dashed hne), and V20R (dash- 
dotted Una) from July 1996 to December 2008. (b) a (solid line) and i?z (dotted line), (c) W (solid line), 
Pq (dotted line), and Pm (dashed line), (d) M (solid line) and E (dotted line). 



Table 1: Correlation coefficients between CME parameters. 



r Vi V20 V2 F20R -a W Po M E R. 



Vi 


1. 


0.998 


0.998 


0.956 


0.670 


0.929 


0.771 


0.724 


0.887 


0.912 


0.894 


V20 


0.998 


1. 


0.994 


0.942 


0.696 


0.932 


0.742 


0.690 


0.893 


0.921 


0.908 


V2 


0.998 


0.994 


1. 


0.967 


0.637 


0.923 


0.794 


0.749 


0.877 


0.899 


0.888 


V20R 


0.956 


0.942 


0.967 


1. 


0.436 


0.820 


0.886 


0.853 


0.768 


0.794 


0.768 


—a 


0.670 


0.696 


0.637 


0.436 


1. 


0.805 


0.205 


0.147 


0.791 


0.793 


0.795 


W 


0.929 


0.932 


0.923 


0.820 


0.805 


1. 


0.630 


0.588 


0.926 


0.892 


0.943 


Po 


0.771 


0.742 


0.794 


0.886 


0.205 


0.630 


1. 


0.994 


0.569 


0.563 


0.551 


Pm 


0.724 


0.690 


0.749 


0.853 


0.147 


0.588 


0.994 


1. 


0.518 


0.502 


0.485 


M 


0.887 


0.893 


0.877 


0.768 


0.791 


0.926 


0.569 


0.518 


1. 


0.971 


0.936 


E 


0.912 


0.921 


0.899 


0.794 


0.793 


0.892 


0.563 


0.502 


0.971 


1. 


0.916 



suggesting that wider CMEs tend to decel- 
erate faster. 

5. Both Pm (dotted line) and Pq (dashed line) 
in Fig. Ic are well correlated with Pi, V20, 
V2, and V20R, r = 0.690 to 0.886, suggesting 
that faster CMEs tend to be closer to the 
solar equator. 

6. Pin is highly correlated with Pq, r = 0.994, 



with the former slightly larger than the lat- 
ter (5.1°). 

7. W is positively correlated with Pq and Pm, 
r = 0.630 and 0.588, respectively, suggesting 
that wider CMEs tend to be closer to the 
solar equator. 

8. a is almost uncorrelated with both Pq and 
Pm, r = -0.205 and -0.147, respectively. 
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9. M (solid line in Fig. Id) is well correlated 
with Vi, V20, V2, and V20R, r = 0.768 to 
0.893, suggesting that faster CMEs carry 
more mass outward. 

10. M is highly correlated with W, r = 0.926, 
suggesting that wider CMEs carry more 
mass outward (Gopalswamy 2010; Aarnio 
et al. 2011). 

11. M is correlated with Pq and Pm, r = 0.569 

and 0.518, respectively, suggesting that 
CMEs closer to the solar equator tend to 
carry more mass outward. 

12. E (dotted line in Fig. Id) is highly correlated 

with M, Vi, V20, V2, and V20R since E = 
(Vourlidas et al. 2000). 

13. E is well correlated with W, r = 0.892, 
suggesting that wider CMEs carry away 
more energy (Gopalswamy 2010; Aarnio et 
al. 2011). 

14. E is correlated with Pq and P,jj, r = 0.563 
and 0.502, respectively, suggesting that 
CMEs closer to the solar equator tend to 
carry away more energy. 

15. The correlation coefficients of Pm arc sys- 
temically lower than those of Pq with other 
parameters, which may be due to a more 
non-linear behavior close to the leading edge. 

16. The parameters arc well correlated with 
(dotted line in Fig. lb), except for an anti- 
correlation of a with Rz. 

17. a is well anti-correlated with M (r = 
-0.791). 

In summary, CMEs with faster speeds tend to 
be wider, to be closer to the solar equator, to de- 
celerate faster, and to carry more mass and energy 
outward. For other properties of CMEs, the read- 
ers can refer to excellent reviews in the literature 
(Gopalswamy 2006; Cremades & St. Cyr 2007; 
Gopalswamy 2010; Gerontidou et al. 2010). 

Aarnio et al. (2011) pointed out that the mass 
of a CME is unrelated to its acceleration. How- 
ever, it is shown that a is well anti-correlated with 
M (r = —0.791). In fact, a is highly correlated 
with the reciprocal of the mass 1/M (r = 0.899), 
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Fig. 2. — (a) a (solid line) is well correlated with 
1/M (dotted line), r = 0.899. (b) Ma (solid line) 
is well anti-correlated with R^ (dotted line), r = 
-0.750. 



as shown in Fig. 2a. It seems to suggest that the 
CME motion approximately obeys Newton's sec- 
ond law. Thus, F = Ma can be understood as 
the 'force' to drive a CME. This force is found to 
be well anti-correlated with R^ (r = —0.750), as 
shown in Fig. 2b for Ma (solid line) and R^ (dot- 
ted line). P > around the solar minima, P = 
around Pz ~ 41, while there is not a definite sign 
of P around the solar maximum. It is well known 
that there are two peaks in R^ for Cycle 23, the 
first one (120.8 in April 2000) being higher than 
the second one (115.5 in November 2001). One 
can see in Fig. 2b that there are also two peaks in 
-(Ma) for Cycle 23, the first one (4.6 x 10^^ in 
December 1999) being lower than the second one 
(5.9 X 10^^ in May 2002). Around the two peaks in 
Rz, F < 0, while coincident with the gap between 
both peaks, P > 0. 
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3.2. Relationships between CME Parame- 
ters and -Rz Described by an Integral 
Response Model 

It is well known that the occurrence of the 

largest CMEs tcmds to peak some years after the 
sunspot maximum (Gonzalez & Tsurutani 1987; 
Tsurutani et al. 2006; Ramesh 2010; Kilcik et al. 
2011). To explain the time delays between CME 
parameters with respect to Rz, in this section, we 
employ the following integral response model (Du 
2011a), hereafter Model I, to analyze the relation- 
ships between some typical CME parameters (V20, 
W, Po and E) and R^, 



y{t) = Dj^,_^x{t')e-(*-''y^dt'+yo, 
(Model I) 



(1) 



where y represents one of the CME parameters 
and a; = iJz, yo is a constant, D is called the 
'dynamic response factor', and r is called the 're- 
sponse time scale' of y to x. Figure 3a illustrates 
the reconstructed series (l^of) dotted line) based 
on this model for y = V20 (solid line) and x = R^ 
(dashed line). 

It is apparent in Fig. 3a that the reconstructed 
series (V20f) reflects well the profile of ¥20- The 
correlation coefficient between V20 and V2of {ri = 
0.966) is higher than the linear correlation coeffi- 
cient between V20 and R.^ {ra = 0.908). The fitted 
parameters of this model are listed in Table 2, in 
which a refers to the standard deviation. 

Figure 3b depicts the reconstructed series {Wf, 
dotted line) based on Model I for y = W (solid 
line) and x = R^ (dashed line). The correla- 
tion coefficient between W and Wf (rf = 0.950) 
is slightly higher than that between W and R^ 
(ro = 0.943). Similarly, the correlation coeffi- 
cient between Pq and the reconstructed series Pgf 
(dotted line) based on Model I (rf = 0.679) is 
higher than that between Pq and R^ (ro = 0.551; 
see Fig. 3c). The correlation coefficient between 
E and the reconstructed series {Ef, dotted line) 
based on Model I (rf = 0.947) is higher than that 
between E and R^ (ro ^ 0.916, see Fig. 3d). The 
above results are listed in Table 2, in which the 
last row indicates the relevant averages. The av- 
erage 'response time scale' is r = 10.8 months. 
The correlation coefficients increased from 0.830 
to 0.886 in average when using Model I. There- 
fore, Model I can better describe the relationships 



between CME parameters and R^. 

In Model I, yo is a constant, representing the 
part of CMEs that are uncorrelated with R^ (re- 
lated to variability of other solar phenomena), D is 
the 'dynamic response factor' of y (a CME param- 
eter) to x{Rz), representing the initial efficiency 
of CMEs {dy/dx\t'=t), and r is the 'response time 
scale' of y to x, representing that an input (so- 
lar activity) may affect the output in the subse- 
quent time period of about r (months) accord- 
ing to an exponential decay factor (e~(*~* ^^^). 
It implies that CMEs are related to the previ- 
ous accumulation of solar magnetic energy. As 
is well known, the upper chromospheric activity 
indices tend to lag behind the sunspot number by 
one to several months, depending on the index, 
which was interpreted in terms of active regions 
evolving from the photosphere upward (Bachmann 
& White 1994; Aschwanden 1994; Temmer et al. 
2003). Various solar magnetic activities evolve 
from the photosphere to upper chromosphere with 
different speeds and times. On average, the long- 
term evolution times (or periodicities) and the 
above lag times are reflected in the 'response time 
scale' (r w 11 months). 

3.3. Relationships between CME Petrame- 
ters and R^ via a Double Decay Inte- 
gral Response Model 

As there are often two (or three) peaks in 
several parameters associated to solar activity 
(Gnevyshev 1967; Feminella & Storini 1997) or 
CMEs (Gonzalez & Tsurutani 1987; Yashiro et al. 
2004; Tsurutani et al. 2006; Kane 2006; Cremades 
& St. Cyr 2007; Ramesh 2010), we use the follow- 
ing double decay model (Model II, hereafter). 



yit) = /;_oo-(^' 

(Model II) 



dt' + yo, 



(2) 



to better describe the relationships between the 
CME parameters and R^. The results based on 
Model II are shown in Fig. 4 and Table 3. 

One can see in Fig. 4a that the correlation co- 
efficient of y = V20 (solid line) with the recon- 
structed series l^of (dotted line) from x ^ R^ 
(dashed line) by Model II is now rf = 0.980, 
slightly higher than that of V20 with the recon- 
structed series based on Model I in Fig. 3a (0.966). 
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Fig. 3.- (a) V20 (solid line), R.^ (clashed line), and the reeonstructed series (V20f, dotted line) by Model 
I (Equation (1)). The correlation coefRcients of V20 with and V20f are ro = 0.908 and rf = 0.966, 
respectively, (b) Similar results for the relationship between W (solid line) and R^ (dashed line). The 
correlation coefficients of W with R^ and the reconstructed series Wi (dotted line) are ro = 0.943 and 
rf = 0.950, respectively, (c) Relationship between Pq (solid line) and R^ (dashed line). The correlation 
coefficients of Pq with R^ and the reconstructed series Pof (dotted line) are ro = 0.551 and rf = 0.679, 
respectively, (d) Relationship between E (solid line) and R^ (dashed line). The correlation coefficients oiE 
with i?z and the reconstructed series Ei (dotted line) are ro = 0.916 and rf = 0.947, respectively. 



Table 2: Fitted parameters of the integral response model (Equation (1)). 



y-x 


D 


r 


yo 


a 


ro 


rf 


V20--Rz 


1.66 


8.8 


549 


159.9 


0.908 


0.966 


W-R, 


0.558 


3.1 


87.0 


26.9 


0.943 


0.950 




0.0381 


24.7 


150 


31.4 


0.551 


0.679 


E/IQ^°-R^ 


0.0133 


6.6 


-0.09 


1.26 


0.916 


0.947 


Average 




10.8 






0.830 


0.886 



The lag time of V20 with respect to Rz {Lq = 28 
months) is well predicted by this model {L{ = 23 
months). 

Figure 4b illustrates the results ioi y = W 
(solid line) and x = Rz (dashed line). The correla- 
tion coefficient of W with the reconstructed series 
Wi (dotted line) by Model II is now rf = 0.954, 
slightly higher than that of W with the recon- 



structed scries based on Model I in Fig. 3b (0.950). 
About half of the lag time of V20 with respect to 
Rz {Lo = 7 months) is predicted by this model 
(Lf = 3 months). In Fig. 4c, the correlation co- 
efficient of Pq (solid line) with the reconstructed 
series Pof (dotted line) from Rz (dashed line) by 
Model II is now rj = 0.734, slightly higher than 
that of Po with the reconstructed series based on 
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Fig. 4. Similar to Fig. 3 but using Model II (Equation (2)). (a) The correlation coefficient of V20 (solid line) 
with the reconstructed series V2of (dotted line) from (dashed line) by Model II is now rf = 0.980. The 
lag time of V20 (^of ) with respect to at the peak of Cycle 23 is Lq = 28 {L{ = 23 months) months, (b) 
The correlation coefficient of W (solid line) with the reconstructed scries Wf (dotted line) from R.^ (dashed 
line) by Model II is now rf = 0.954. The lag time of W (Wf ) with respect to R^ is Lq = 7 {Lf = 3 months) 
months, (c) The correlation coefficient of Pq (solid line) with the reconstructed series Pof (dotted line) from 
Rz (dashed line) by Model II is now rf — 0.734. The lag time of Pq (Pof) with respect to R^ is Lq ~ 25 
{L{ = 27 months) months, (d) The correlation coefficient of E (solid line) with the reconstructed series E{ 
(dotted line) from R^ (dashed line) by Model II is now rf = 0.954. The lag time of E (Ef) with respect to 
i?z is Lq = 26 {Lf = 22 months) months. 



Table 3: Fitted parameters of the double decay integral response model (Equation (2)). 



y-x 


Di 


ri 


D2 


T2 


yo 


(T 




r-f 


Lo 


Lf 


V20-R. 


6.30 


1.6 


0.13 


75.6 


286 


123.5 


0.908 


0.980 


28 


23 


W-R. 


0.94 


1.4 


0.0092 


34.9 


74.9 


26.6 


0.943 


0.954 


7 


3 


Pq-Rz 


0.0811 


5.0 


0.0167 


82.9 


121 


29.8 


0.551 


0.734 


25 


27 


E/10^°-R^ 


0.0092 


17.0 


-0.0022 


45.3 


1.12 


1.17 


0.916 


0.954 


26 


22 


Average 




10.8 






0.830 


0.886 











Model I in Fig. 3c (0.679). The lag time of 1^20 

with respect to R.^ (Lq — 25 months) is approxi- 
mately predicted by this model {Lf = 27 months). 
Similarly in Fig. 4d, the correlation coefficient ofE 
(solid line) with the reconstructed series Ef (dot- 
ted line) from R^ (dashed line) by Model II is now 



Vf = 0.954, slightly higher than that of E with the 
reconstructed scries based on Afodcl I in Fig. 3d 
(0.947). The lag time of V20 with respect to i?z 
{Lq = 26 months) is approximately predicted by 
this model {Lf = 22 months). 

In Table 3, the last row indicates the relevant 
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averages. The average correlation cocfScicnt based 
on Model II (ff2 = 0.906) is slightly higher than 
that based on Model I (ffi = 0.886) and higher 
than the linear correlation coefficient {Tq = 0.830). 
The two average response time scales of Model II 
are about ri = 6.3 and T2 = 59.7 months, one 
being shorter and the other one being longer than 
that for the case of Model I (r = 10.8 months). 
Model I is therefore a simplified version of Model 
II. 

It has been known that low-energy phenom- 
ena {e.g., sunspots, 10.7-cm radio flux and Ca 
K index) are associated with lower atmospheric 
layers {e.g., photosphere, chromosphere and up- 
per chromosphere, respectively), while high en- 
ergy phenomena {e.g., CMEs) are associated with 
higher atmospheric layers {e.g., corona) and origi- 
nate mostly from the solar active regions (Gosling 
et al. 1976; Bachmann & White 1994; Feminella & 
Storini 1997; Sheeley et al. 1999). The solar mag- 
netic structures (represented by some activity in- 
dices) may be grouped into two classes: short-lived 
and long-lived. Similarly to the case described in 
the previous section, CMEs are assumed to be re- 
lated to the previous accumulation of magnetic en- 
ergy, but now with two 'response time scales': the 
first (ti « 6 months) is related to short-lived weak 
structures which tend to lag with respect to i?z by 
several months (Bachmann & White 1994; Tem- 
mer et al. 2003), and the second {t2 « 60 months) 
is related to long-lived strong ones which tend to 
peak later than by about a few years (Aschwan- 
den 1994; Bromund et al. 1995) and may be related 
to the 5.5-year periodicity in solar activity (Das & 
Nag 2003). 

4. Discussions and Conclusions 

It has been shown that the CME parameters 
(speed, acceleration, central position angle, angu- 
lar width, mass and kinetic energy) are well corre- 
lated, with high correlation coefficients (from 0.77 
to 0.99) between most of them. These indicate 
that CMEs with faster speeds tend to be wider, 
to be closer to the solar equator, to decelerate 
faster, and to carry more mass and energy out- 
ward. After 2004, V20R becomes larger than the 
other three speeds (Vi, V20 and V2), which is re- 
lated to the increasing trend in the acceleration a 
(solid line in Fig. lb). In addition, a is found to 



be highly correlated with the reciprocal of mass, 
1/M {r = 0.899), suggesting that the CME motion 
obeys Newton's second law. The force to drive a 
CME is found to be well anti-correlated with 
{r = -0.750). 

The relationships between some typical CME 
parameters (V20J Pa and E) and sunspot num- 
bers (i?z;) can be well described by an integral 
response model (Equation (1)) — Model I. The 
correlation coefficients between the CME param- 
eters and Rz increase by 6.7% from To = 0.830 
to ffi = 0.886 on average when using this model 
(Table 2). In this model, the output y{t) depends 
not only on the present input x{t) but also on 
its past values according to an exponential decay 
factor The earlier the input, the less 

it contributes to the output. This implies that 
CMEs are related to the previous accumulation 
of solar magnetic energy with a mean 'response 
time scale' of about r « 11 months. Apart from 
short-term variations, parameters related to solar 
magnetic activity have also long-term components 
evolving from the photosphere to the upper chro- 
mosphere. The magnetic energy (in some events) 
can be stored to be released later instead of being 
released instantaneously. As the CME parameters 
are well cross-correlated, similar conclusions can 
also be drawn when using the other six CME pa- 
rameters in Table 1. 

Using the double decay integral response model 
(Equation (2)) — Model II, the above correla- 
tions improve (rf2 = 0.906). These results mean 
that about Tq — 68.9% (coefficient of determi- 
nation) of the variations in the CME parameters 
can be explained by the variation in Rz through 
a linear relationship. While about ~ 78.5% 
{7(2 = 82.1%) of the variations in the CME pa- 
rameters can be explained by the variation in R.^ 
through Model I (II). Besides, the time delays of 
CME parameters with respect to R^ {Lq = 22 
months) can also be well predicted by the model 
(Lf = 19 months), Lf/Lo = 86%. Therefore, 
Model I or II can better describe the relationships 
between CME parameters and Rz. 

Model II implies that the relationships between 
CME parameters and sunspot (magnetic field) ac- 
tivity depend on two exponential decay factors, 
one with a shorter time scale (ti = 6.3 months) 
and another one with a longer time scale (t2 = 
59.7 months). This model is related to the two (or 
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three) peaks often present in solar activity indi- 
cators (i?z, flares, radio, and X-ray fluxes, Gnevy- 
shev 1967; FemineUa & Storini 1997), CMEs (Gon- 
zalez & Tsurutani 1987; Tsurutani ct al. 2006; 
Ramesh 2010; Kilcik et al. 2011), geomagnetic ac- 
tivities or CME-related storms (Gonzalez & Tsu- 
rutani 1987; Tsurutani et al. 2006), one near the 
peak in and another one a few years later. A 
double peak in these activity indicators suggests 
that there are two sources or two decay time scales. 

It has long been recognized that low-energy 
phenomena (slow CMEs) are associated with 
deeper atmospheric layers (erupting prominences) 
and tend to follow the sunspot activity, while high 
energy phenomena (fast CMEs) originate mostly 
from active regions and tend to lag behind the 
sunspot activity (Gosling et al. 1976; Feminella & 
Storini 1997; Sheeley et al. 1999). The two 're- 
sponse time scales' in Model II (ti w 6 months 
and T2 « 60 months) may be related to the short- 
lived weak magnetic structures (low energy activ- 
ity phenomena) which tend to lag behind by 
several months (Bachmann & White 1994; Tem- 
mer et al. 2003) and the long-lived strong ones 
which tend to peak later than i?z by about 2-3 
years (Aschwanden 1994; Bromund et al. 1995). 
This reflects the fact that solar magnetic activity 
evolves from the photosphcire to the upper chro- 
mosphere with different speeds and different time 
scales (or periodicities). The 'response time scale' 
(r « 11 months) in Model I is the mean effect of 
the two (ti w 6 months and T2 « 60 months) in 
Model II. 

Solar structures near the solar equator are more 
active than those near the solar poles, which may 
be related to the faster solar rotation at lower 
latitudes. According to the above model(s), the 
stronger the inputs (x), the more they contribute 
to the outputs (y), and the longer the lag time of y 
with respect to x. Large active regions of complex 
magnetic structure are long-lived and produce in- 
tense CMEs which have long time delays with re- 
spect to Rz. Therefore, active regions close to the 
solar equator tend to generate CMEs with faster 
speeds which carry more energy outward. 

As pointed out in Section 2, the above results 
may be affected by CME projection effects, be- 
cause the CME parameters are measured in the 
plane of the sky. More accurate results could be 
obtained using data from the Solar Terrestrial Re- 



lations Observatory (STEREO) as we can deduce 
the three-dimensional nature of CMEs. As full 
halo CMEs account only for 3.6% of all CMEs 
and CMEs with W > 120° accoimt for 11% of 
all CMEs (Gopalswamy et al. 2010), and in the 
present study we use the integrated values of CME 
parameters, we expect that our general conclu- 
sions are not significantly influenced by projection 
effects. 

The main points of this study can be summa- 
rized as follows, 

1. The acceleration of a CME (a) is highly cor- 
related with the reciprocal of its mass (M), 
r = 0.899. 

2. The force (Ma) to drive a CME is well anti- 
correlated with the sunspot number (i?z), 
r = -0.750. 

3. The relationships between CMEs and R^ can 
be well described by an integral response 
model (Equation (1)) with a decay time scale 
of about 11 months. The correlation coeffi- 
cients of CME parameters with the recon- 
structed series based on this model (ffi = 
0.886) are higher than those based only on 
i?z {ro = 0.830). 

4. A double decay integral response model 
(Equation (2)) with two decay time scales 
of about 6 and 60 months shows better cor- 
relations between CME parameters and Rz 
(ff2 = 0.906). 

5. The time delays of CME parameters with 
respect to Rz can also be well predicted by 
the model (86%). 
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